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Abstract 

In a recent paper we derived the free energy or partition function of the A^-state chiral Potts 
model by using the infinite lattice "inversion relation" method, together with a non-obvious 
extra symmetry. This gave us three recursion relations for the partition function per site Tpq 
of the infinite lattice. Here we use these recursion relations to obtain the full Riemann surface 
of Tpq. In terms of the tp,tq variables, it consists of an infinite number of Riemann sheets, 
each sheet corresponding to a point on a (2A^ — l)-dimensional lattice (for N > 2). 

The function Tpq is meromorphic on this surface: we obtain the orders of all the zeros 
and poles. For TV odd, we show that these orders are determined by the usual inversion and 
rotation relations (without the extra symmetry), together with a simple linearity ansatz. For 
A'^ even, this method does not give the orders uniquely, but leaves only [(A + 4)/4] parameters 
to be determined. 

Keywords: Statistical mechanics, lattice models, chiral Potts model, free energy 

1 Introduction 

The free energy of the A-state chiral Potts model was obtained implicitly in 1988 from the 
star-triangle relation, and explicitly in 1990 [^], Q from the functional relations for the finite- 
lattice transfer matrices. Recently |^ we considered the infinite-lattice limit of these functional 
relations, and showed that they reduced to the rotation and inversion relations (which exist 
for most two-dimensional models, solvable or not |^), or rather to a single combined rotation- 
inversion relation. All other equations were either definitions of the auxiliary function T2{tq), 
or direct consequences of the relations and definitions, so by themselves contained no new 
information. 

We showed that the reason we could solve the infinite relations was that we knew the 
analyticity properties in some "central" domain that included part of the physical regime. In 
particular, we knew that T2{tq) had only one branch cut in the tg-plane, instead of the N cuts 
one might expect. This was equivalent to an extra symmetry for the free energy. It ensured 
that one could use Wiener- Hopf methods to calculate the partition function per site Tpq. 

We found that Tpq satisfied three "recursion" relations; equations (70), (71) and (75) of 
(We could have obtained these directly from the solution found in 1990.) Here we use these 
three relations to obtain the full Riemann surface (in both the p and the q rapidity variables) 
on which Tpq lives. It is a meromorphic function on this surface, its only singularities being 



±1 

<3 



poles. Its poles and zeros occur only when tp — tq , Xp — X, 

Regarding Tpq as a function of tp 
of Riemann sheets. Each sheet is specified by a set of N integers mo, . . . ,mM-i associated 
with the p-variable, and another set no,...,n]v-i associated with the g-variable. Thus it 
corresponds to a point on a 2A'^-dimensional lattice. In fact, although the chiral Potts model 
does not have the usual difference property (that one can choose the rapidites p, q so that the 
Boltzmann weights, correlations and free energy depend only on p, q via their difference p — q), 
we do find a weak residual version of this property: Tpq is the same on all sheets obtained by 
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incrementing mo, . . . , mjv-i, ^^0l • • • , n^-i by an arbitrary integer. Hence the Riemann surface 
can be associated with a lattice of reduced dimension 2N — 1. 

The case N = 2, when the model reduces to the Ising model, is special. Then two of 
the integers mo, . • . ,njv-i vanish, so the surface is one-dimensional, being that of the elliptic 
function argument z = exp[n{up — Uq)/K'] of equation (B7) of 

A knowledge of the Riemann surface, and of the orders of the poles and zeros, gives us a 
great deal of information about the function Tpq. We observe that that these orders are linear 
in the mo, . . . ,mM-i, and in the no, ... , um-i- Hence they are bilinear in the full set of 2N 
integers mo, . . . ,mjv-i,no, . . . ,nN-i, depending only on their differences. 

We ask the question: is this bilinearity, together with the basic inversion and rotation 
relations sufficient to determine the orders of the poles and zeros? We find that for A'' odd 
the answer is yes. For N even it is not quite sufficient: there are [{N + 4)/4] parameters still 
undetermined. 



2 Domains and Riemann sheets 

The Boltzmann weights of the A-state chiral Potts model are functions of a constant k, and 
two "rapidity" variables p and ^ Let 

Then the variable p can be thought of as a point {xp, yp,tp,\p) on the algebraic curve 

■^p ' Vp — "'V-'- ' -^p Vp ) : '^p — ^pUp 5 

k = 1 — k' /\p , ky^ = 1 — k' \p . (2) 

Similarly q is the point {xq,yq,tq). 

As we remark in if Xp,Xq,yp,yq,L0Xp all lie on the unit circle and are ordered anti- 
cyclically around it, then all the Boltzmann weights are real and positive, and therefore so is 
the partition function per site Tpq. We refer to this case as the physical regime. Outside this 
regime we define Tpq by analytic continuation. 

This function Tpq therefore lives on a Riemann surface in both the p and the q variables. 
To specify this surface, first consider the p variable. If |Ap| < 1, then Xp lies in the region 
S of Figure 1, while j/p lies in one of the N approximately circular regions TZo, . . . ,TZn-i 
surrounding the points . . . ,u!^~^. 

These regions shrink to points in the low-temperature limit k' 0, so in this limit it is 
certainly true that 

yp ~ Lo'^ if yp e TZr ■ 

We find to helpful to write " yp ~ ci;'" " , by which we mean " yp £ TZr " . 

We define a "domain" Dr for Xp and yp simultaneously by saying that if yp lies in TZr , 
then {xp, yp) lies in the domain T?r. More simply, we say that p lies in T>r. We also say in this 
case that p has parity and that it has type r. 

Conversely, if |Ap| > 1, then j/p lies in the region S and Xp lies in one of TZo, . . . ,TZn-i- If 
Xp lies in TZr , then we say that {xp, yp) or p lies in the domain T>'r, has parity 1 and type r. 

Let p, r be the parity and type of p. Then 

Xp € S , yp eTZr , |Ap| < 1 and p e Vr if p = , 

Xp £TZr , yp € S , |Ap| > 1 and p G if p = 1 . (3) 

We refer to the case when p and q both lie in Do, so that |Ap| < 1, lA^I < 1, yp,yq £ TZo, 
the central regime or domain. It overlaps the physical regime, so Tpq is readily 
extended to this domain. Series expansions strongly suggest that it has no poles or zeros in 
the central domain (as can be verified from the explicit solution (67) of Q). 

If we analytically continue from the central domain to |Ap| > 1, then as Ap crosses the unit 
circle, Xp enters one of the regions TZo, ■ ■ ■ ,TZn-i, say TZn, while j/p leaves TZo and enters S. 
Thus p goes from the domain I?o to D'r-^ . Its type changes from to ri, and its parity from 
to 1. 
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We can then analytically continue from |Ap| > 1 back to |Ap| < 1, but p will not necessarily 
return to Cq: in general it will go to a new domain ©rj- And so on: p will move successively 
through a sequence of domains 

©0 , -D'r, , T^r, , , I'ri , • ■ ■ (4) 

Their types are 0, ri, r2, ra, r4, . . ., and their parities are 0, 1, 0, 1,0,.. .. 
From (|), 

k^t^ = l-fc'(Ap + l/Ap) + fc'^ , 

so the unit circle in the Ap-plane corresponds to N straight-line segments in the tp-plane, from 
ut^rj to /ri, as in Figure 2 of [^. Here 77 = [(1 — fc')/(l + k')]^^'^ . In terms of the tp variable, 
each {xp,yp) domain is therefore a Riemann sheet, with branch cuts along these straight-line 
segments. As one goes from a domain to a neighbouring domain, tp crosses one of these A'' 
branch cuts and goes from one sheet to the next. We shall use the words "domain", "Riemann 
sheet" and "sheet" interchangeably. 

We define Tpq by analytic continuation as p moves through this sequence, starting from 
some initial value in Do- As p moves from one domain to the next, Tpq moves from one 
Riemann sheet to the next. 

At this stage we do not know the Riemann surface on which Tpq lives (i.e. how the Riemann 
sheets connect with one another), so we must be prepared to remember this full sequence of 
domains in order to uniquely specify the value of Tpq at a point p. Since the parities alternate, 
it is sufficient to remember the domain types {ri, r2, ra, 7-4, . . .}. 

Similarly, if initially the other rapidity variable q is in the central domain ©o, and is 
analytically continued through the sequence 

1^0 , ^^si , 'C>S2 , T^i^ , Vs^ , . . . 

we must remember the sequence {si, S2, S3, S4, . . .}. 

Fortunately it seems that Tpq does not depend on whether p or q moves first: i.e. how the 
two sequences {ri, r2, fa, r4, . . .}, {si, S2, S3, S4, . . .} are interleaved. Thus the Riemann surface 
in both variables is simply the union of the two single- variable surfaces. 





Figure 1: The iV + 1 regions S, TZq, . . . , TZn-i of the complex plane 
in which Xq and yq lie (for iV = 3). 7?.o, . . . TZn-i are the 
interiors of the approximate circles centred on 1, . . . ,uj^~ 
S is the the complex plane outside all N circles. 
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3 Recursion relations for Tpq 

In we obtained the two relations, true for r, s = 0, . . . , A'^ — 1: 
T (oj^-ii r \t - NT{uj'-Xp,yp\xg,y,) ^ tp - lj-H, 



T{u ^Xp,yp\xg,yq)T{xp,yp\xg,yg) i-i-{xp-uj 3Xq){yp-L0 iyq) 



JV-l 



TT '^i'"^ (5) 

J-i {xp- u)-Oyq){yp- u-^Xq) 

j = r+l 



NT{xp,yp\u"Xq,yq) -Ar tp - 4, 



Tix v\lo'v x) = 1^ iK^p,yp\uj Xq,yq) -pr 

sK v.Vv\ y,, ,) T{xp,yp\u;-^Xq,yq)T[xp,yp\xq,yq)^y[xp-u;^-^yq)[yp-uo^Xq) 

J- J- (Xp - ^Xq){yp-UJ^yq) 

J=S + 1 

together with 

T(j:p,i/p|a::g,i/5)r(a;a;p,yp|cjx,,i/q) {xp ~ Xq){ujtp - tq) 



T{Xp,yp\ujXq,yq) T(ujXp,yp\Xq,yq) (ujXp - Xq)(tp - tq) 



(7) 



In these three relations the p-arguments {u^Xp,yp) (for any i) of the function T (without 
a suffix) lie in the initial domain T>q, and so do the g-arguments {Lu^Xq,yq). The function Tr 
on the Ihs of (^) has p-arguments (x'p,y'p) = (ui^yp, Xp) lying in the domain "D^ adjacent to Do- 
Similarly, Ts on the Ihs of ^ has g-arguments (x'q,y'q) = {uo^yq^Xq) lying in the domain D^ 
adjacent to Do- 

We can use (^ and (^ to obtain the function T on any of its infinitely many Riemann 
sheeets. For instance, if in ^ we allow p = {xp,yp) to leave and enter the neighbouring 
domain I?^, where Xp ~ a;™, then the arguments of T on the Ihs will move to a sheet T>rm- 
The three Ts on the rhs move to the domains T>'m+ri I^m-ij T^'m, respectively. These can in 
turn be expressed in terms of T in Vq by again using 

We can repeat this procedure ad infinitem. At each stage the function Tpq on the left-hand 
side of (P) moves to a new domain or sheet, while the ones on the rhs move to domains or 
sheets that have already been expressed in terms of the values of Tpq in the inital domain Vo- 
Similarly, we can use (p|) to obtain Tpq on successive sheets in the g-variable. 

At first sight we will obtain an exponentially infinite Cayley tree, each (p, q) surface having 
2A'" neighbours (A^ neighbours in the p variable, TV in the g- variable). The number of mth 
neighbours of the initial domain (where p,q £ T>q ) will be 2A'' x (2A'' — 1)™"^. 

However, it is not as as bad as that. Consider the function Tpq on a Riemann sheet of 
p-parity p and p-type r, and of g-parity a and g-type s. Thus p, a = 0, 1 and r, s — 0, . . . , N—1. 
Define an associated function A^^ {ap,bp\aq,bq) by 

T{xp,yp\xq,yq) = A°1{xp,uj~''yp\xq,uj~''yq) if p = 0, a = , 

= N/A°l{xp,uj~''yp\yq,uj~''xq) if p = 0, cr = 1 , 

= N/Al°{yp,u}~''xp\xq,uj~''yq) if p = l,cr = , (8) 

= All(yp,uj~''xp\yq,uj~''xq) if p=l,a = l . 

In each case the p-arguments ap,bp of Af^'^{ap,bp\aq,bq) lie in "Do, and so do the g-arguments 
aq,bq. 

Then no matter how many times we iterate (H) and (H), the resulting function Tpq will 
have an associated function of the form 

JV-l 

Arl{ap,bp\aq,bq) = ^^^[{ap ~ lo^ aq)"' {ap ~ Lo^bq)"' (9) 

3=0 

JV-l JV-l 

x{bp-uj^aqf^{bp~ijj^bqY'i{apbp-uj^aqbqy^]x J]^ '^T{uj^ap,bp\J aq^bq)"'^ , 



i=0 j = 
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where the exponents Qj, . . . ,pij are integers and the functions T on the rhs all have p and q 
arguments lying in the initial sheet I?o. 

Further, we can always use (^) to eliminate the T{uj'^ap,bp\ijj-' aq,hq) with i and j greater 
than zero, so we can require that 

Pij =0 if i and j > , (10) 

leaving only poo, pio and po-, > 0) as possibly non-zero. There are then at most SA'^ + 
2A'" — 1 = 7N — 1 non-zero exponents. Equating any two Riemann surfaces with the same 
values of the function Tpq, it follows that any Riemann surface (in both the p and q variables) 
can be specified by at most 7A'^ — 1 arbitrary integers (together with the two- valued integers 
p, (J and the A''- valued integers r, s), so can be regarded as corresponding to a point on a 
(TAf— l)-dimensional lattice. 

In fact the lattice is of lower dimension still, since there are many relations between the 
exponents. For any rational number x, let [x] denote its integer part, so that 

[x] <x <[x] + l , 

and define 

Fi,j = [^] - [^] - [^] . (11) 

Then Fij has value or 1, and is periodic in i,j of period A'^, i.e. Fi+N,j ~ Fi^j+M = Fij. 

The relations between the exponents are conveniently expressed by associating with each 
Riemann sheet two more sets of integers]^ 

m = {mo, . . . ,mjv_i} , n = {no, . . . , njv-i} . 

The set m is varied when the p-domain is changed, and n when the g-domain is changed. If 
m is the set on a sheet of parities p, o and types r, s, and m' is the set m on a p- neighbouring 
sheet of parities p' ,0 and types r',s (the q variables being the same on each sheet), then 
p' — 1 — p and 

m'j = rrir+r' + "Ir + r'-l-l " "^i + P(l — Sr + r',o) ~ Fj,r+r' , (12) 

for j — 0, . . . , N ~ 1, using the periodic convention 

TTij+jv = nij , Vj . (13) 

Similarly, the set n on a sheet of parities p, a and types r, s is related to the set n' on a 
q-neighbouring sheet of parities p, a' and types r, s' by a' = 1 — ct and 

n'j = n^+s' + ris+s'+i - rij + cr(l - 5s+s',o) ~ Fj^s+s' ■ (14) 

These recursion relations, together with the conditions 

mo = • ■ ■ mjv-i — no ■ ■ ■ = njv-i =0 in the initial sheet, where p,q £ I?o 

define the integers mo, . . . , mjv-i, no, . . . , njv-i on all sheets. In particular, it follows that 

JV-1 JV-1 

nij = Nfi~ r , Uj ^ Nu - s , (15) 



j=0 j=0 



p and u being integers. 
Define 



JV-1 

9k ^ ^2 '^i + k'TI'i ■ (16) 



^We originally found the rij by keeping the p variables in ©o and looking at the orders of the zeros of Tpq, 
considered as a function of q. Then rrij = and the rhs of (Ea) becomes 2rt_fe + toj.. In the normalization where 
the Boltzmann weights are free of poles we lose the w^. and the zeros are simply 2{—l)'^n_f^. Similarly, one can 
obtain the m, . 
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Numerical Fortran experiments (in integer arithmetic) for small TV and sheets close to T> 
strongly suggest that 

Oj — fj. — u — rajj^r — rrij+r+i + ris-j-i + ris-j + Stj , 
a'j — —fi — u + mj+r + riij+r+i + a'j , 

(3j = fl + l^- Us-j - Us-j + l + Pj , 

P'j = -ii + u + p'^ , (17) 
7i ~ 2 (7j4-r-s+i ~ 2 (jj+r-s-i — 2 mj+r-s+CT + 2 rzs_r-j+p + 

rrij+r + rrij+r+i — ris-j-i — Us-j + , 
pij — {nir-i+i — mr-i-i)Sj,Q + (us-j+i — ns-j-i)Sifi + pij , 



where 



ii] + [^^^] - - + (1 - 



LivJ ' L JV J L iV J L AT 

7. = [^] + [^]-[i]-[^]-(l-p)a + p(l-a)(-l + 2 5,,,_0 , 

From (|l5|), incrementing r (s) by A'^ increments /i (i^) by 1, so the right-hand sides of ( [l7| ) are 

each periodic in r, s and j, of period A^, as they should be. 

Define 

0P<T = p + o — 2pa , (19) 
so that (j)pcr = if p = o", and ^po- = 1 if p 7^ cr. Our numerical experiments also suggest that, 

= ^/3/ = (iV - 1)(1 - p)a + ^(n,-m,) , 
i i j 

a/ = ^ = (A^ - l)p(l - a) + ^(m, - n,) , (20) 

J i i 

the sums being from j = to j = A'^ — 1. It follows that 

J2{a, + q/ + 7,) = ^(/3, + /?/ + 7,) = , (21) 



^(a,+ft+7.) = + + 



7. J 



We also find that 



^j(a,+ + ft' +7.) = , modiV (22) 



if N is odd or if p = a; if TV is even and p ^ cr, then the rhs of ( p^ ) equals TV/2, modulo TV. 

These relations ( po| ) - (^2|) ensure that no additional external constant factors occur in (^). 
For instance, multiplying Up by will introduce a factor where L = X]j(aj +aj' + 7j). 
From (^), this factor is unity. Similarly for hp,aq, hq. Also, interchanging all p variables with 
the corresponding q variables will introduce a factor ( — 1)^^;'', where 

I = Y{'^]+a/ + l3,+fi/ + 'yj) , J = ^j(aj+a/ + ft +/3/ + 7j) ■ 
i j 

From ( po[ ) and (|2^), this factor is also unity. 

If we ignore the requirement (|l^) ,0, the relations (^8|) are unchanged by the substitutions 

p, a, r, s,p, i^,mj, nj , Qj , , ft , ft , 7^, p^j — » 1 - a, 1 — p, 1 — s, 1 - r, — i^, — ni-j, -mi_j, 
Oj , ft+i , Oj-i , P'j , 7j , 



Presumably there is some variant of this symmetry that preserves the equations (|17[) completely. 
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"Dimension" of the Riemann surface 



We see that the 2N integers mo, . . . , mN-i, no, . . . , njv-i are sufficient to specify the function 
Tpq on any sheet within its Riemann surface, so any sheet can be associated with a point 
in a 2A'^-dimensionaI space. Further, incrementing each of mo, . . . ,njv_i (and therefore also 
v) by unity (or any integer) leaves (^) unchanged, so the space can be reduced to one of 
dimension 2A'' — 1. Ths appears to be a partial analogue for the chiral Potts model of the 
rapidity "diflerence property" that plays such an important role in the simpler models. 

For N > 2 this appears to be the best one can do - each sheet of the Riemann surface cor- 
responds to a point in a (27V — l)-dimensional space. Note however that sheets are neighbours 
if their m-integers satisfy (^, or if their n-integers satisfy (|l^). Hence neighbouring sheets 
do not necessarily correspond to neighbouring points in mo, . . . , n]v-i-space. 

For N = 2 we have the additional relations rrii-p = ni_o- = for all sheets, so mo, mi, no, 
ni enter ( p^ only via nip — ria and the space is one-dimensional, as we observed in ^ . 

For N > 2we note that the exponents Oj, . . . , /3'j, pij are linear in the integers mo, . . . , njv-i , 
but the 7j are linear only separately in mo, . . . , mjv_i and no, ... , njv_i . They are bilinear in 
the full set of 2N integers, due to the occurrence of the Qj, as defined by (^). 

The values of p = {xp, j/p, p, r, mo, . . . , mjv-i} completely specify the point p, not only on 
the algebraic curve (|^), but also on the Riemann surface of Tpq. Similarly, q = {xq,yq,a, s, 
no, . . . ,n]v-i} completely specifies q. We shall refer to the corresponding Riemann sheet as 
the "sheet (m, n) " . 

4 Zeros and poles of Tpq 

For ap,bp,aq,bq in the central domain V, the functions T(uj^ap,bp\uj^ aq,bq) on the rhs of ^ 
are non-zero and analytic. Hence the rhs of (^ is meromorphic in V, with zeros or poles only 
when 

Qp = LJ-' Uq and bp — bq , (23) 
which implies apbp — uj-'aqbq. Thus its zero at this point is of order 

aj+P'o+lj ■ (24) 

(Equivalently, its pole is of order —aj — P'q — 7j.) 

From (^, the relation Xp — x^ implies i/p — and vice- versa, and either implies 
tp = tq . Similarly, Xp = yq implies yp = x^ and again tp = tq . From (^ it follows that 
the function T{xp,yp\xq,yq) is meromorphic throughout its Riemann surface, with zeros and 
poles only when tp = . On the sheet (p, a, r, s, m, n) the zero at tp — uj'^tq (fc = 0, . . . , A''—!) 
is of order 

Zk = { — ly*" {OLk + s^r + P'o + "fk + s-r) ■ 

From (|l^ and (|l|) it follows that 

{-lY+''zk = 2(7fe+i-2gfc_i-2mfc+<,-f 2np_fc+2p(l-a)(5fe,o-l) + Wfc , (25) 

where 

Wk = [ "+-'■-^1--) ] _ [^] + + . (26) 

The contribution Wk has a simple explanation. Using the notation of [^] and [^, let us 
postulate the existence of functions Oi,j, Qi.j of p and q such that Qi^j has simple zeros on 
the Riemann surface when 

Xq = uj^yp and yq = uo'' Xp , (27) 

and 0ij has simple zeros when 

Xq = u^Xp and yq = u-'yp . (28) 

Any zero or pole of Tpq must occur at one of these points on some particular Riemann sheet. 
Define Gpq by 

Af-l JV-i 

Gpq = G{Xq,yq) = H ^ ,N +1- jQ ^ ' , (29) 
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Then Gpq is the function 1/^^ of We remark therein that its poles contain just all the 
poles of the Boltzmann weight functions Wpq{n),Wpq{n). Hence for a finite square lattice of 
A/" sites with partition function Z, it is true that Z/G^g has no poles on the Riemann surface. 
It is therefore natural to define a normalized function Tpq by 

Tpq — GpqTpq , ^30) 

since this Tpq is the partition function per site in this pole-free normalization. 

On the sheet {p,a,r,s,m,n) we find from ( ps| ) that the pole of Gpq at tp = ui'^tq {k = 
0, . 1) is of order —( — 1)'''*''^ Wk (so numerically this is either or 1). Hence the term Wk 

in is the contribution of Gpq. The preceding terms (all even integers) are the contribution 

of Tpq . 

Let 9k be the function with simple zeros on the Riemann surface when tp — Lu'^tq and 
< = <, < = . Then 

iV-l 

ffc ~ n e,,_fc_, , (31) 

i=0 

where by / ~ <; we mean herein that / and g have the same zeros and poles (when tp = tq) 
on the Riemann surface. 
Similarly, 



H'^'.-fc-' ' (32) 



where Ok has simple zeros when tp = uj''tq and Xp = Xq, yp — . Then formally we can 
write (^) as 



Wiek/ekY'"^''' , (33) 



where the product is over all Riemann sheets (p, iJ,r,s, m, n) and all values 0, . . . , A'^ — 1 of fc, 
and the integer tpq{k) is 

€pq{k) = gk+1 - gk-1 - ruk+cr + rip-k + p{l - cr){5k,o - 1) ■ (34) 

If p = (T then only 9k can have zeros, and if p 7^ a only 9k- 

We note that there are various formal identities between our O and functions, notably 



6ij ~ j/g - oj^^p , ]^6ife6lfe ~ tp - tj'^f, , (35) 

i k 

all products being over the integers 0, . . . , A'' — 1. 

It may be possible to give explicit representations of our postulated functions Qij , Qij , Gpq , 
Ok, 9k, Tpq in terms of hyperelliptic functions |[ ^], but we shall not do so here. Using them 
does greatly simplify the relations and certainly provides a way of keeping track of the poles 
and zeros at tq — tp on the Riemann surface. These are the only poles and zeros, apart 
possibly from ones occuring when p has some particular value independent of q, or q has some 
value independent of p. This suggests that herein the relation 

fvi ~ dvi 

is equivalent to the explicit identity 

fpqfp'q' _ 9pi9p'q' 



fpq'fp'q 9pq'9p'q 

for all rapidities p,q,p',q'. This implies that there exist single-rapidity functions Up,Vq such 
that 

fpq = UpgpqVq . 

In all the cases where we have been able to test this hypothesis, e.g. by using (^^, we have 
found it to be true. 
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5 Automorphisms 

Three basic automorphisms that take a point on on the curve to another such point are 
R, S, U, where 

R : XRp = yp , yRp = cuxp ; S : xsp = yp^ , ysp = x~^ ; (36) 

U : xup — ioxp , yup = yp ■ 

lip hes in the initial domain (so ~ 1), then we can take Up to also lie in T>(,\ and Rp, 
Sp to lie in the adjacent domain ©oPI 

We can determine what happens when p is not in the initial domain Do by analytic con- 
tinuation. If p lies on a sheet that is a kih neighbour of Do, a route to it being the sequence 
(^) of sheets of types {0,ri,r2, . . . ,rk}, then Up is also on a fcth neighbouring sheet (so has 
the same parity p), but from ( |36| ) the sequence to Up is {0, ri + 1, r2, rs + 1, . . . , + p}|^ 

As above, writing p = {xp,yp, p,r,mj} for the full set of parameters defining the rapidity 
p (with J = 0, . . . , A'' — 1), it follows that 

Up = {ujXp,yp,p,r + p,mj-i + {-l)''{mN-i - mo) - pSj^i} . (37) 

Iterating , we obtain for all integers i 

U^p = {ijj'xp,yp, p,r + ip,mj-i + {-lYiniM-i - mo) - pFj^i} . (38) 

The r values should always be taken modulo A'^, so that < r < A*'. Then we see, using (|l^), 
that U'^p = p. 

Similarly, Rp and Sp lie on {k + l)-th neighbouring sheets of Do, at the termini of the 
routes {0, 0, n + 1, r2, rs + 1, . . . , + p}, {0, 0, -n, -r-z, -rz, . . . , -r^}, respectively, and 

Rp = {Vv^^v^ - P^r + p,mj-i+ p{l-5j,x)} , (39) 

Sp = {l/yp,l/xp,l ~ p,-r,-mN+i-j} ■ (40) 
Two combinations of automorphisms that we shall need are V — RU^^ and U^V: 

Vp = RU'^p = {yp,Xp,l - p,r,mj + {-l)''{mx - mo)} (41) 

and 

C/'Vp = {Lo''yp,Xp,l - p,r + i{l ~ p),mj-i + {-l)''{mi ~ mN-i) ~ {I - p)Fj,i} (42) 
for any integer i. 

The efi'ect of these automorphisms on the second rapidity q can be obtained at once by 
replacing p, p,r,mo, . . . ,mjv_i by q,a, s,no, . . . ,nN-i- 

6 Relations for Tpq and its exponents 

Going to the pole- free normalization (^o[), the relations (^, (^), (^ simplify to 

rf, , r I X T{uj''Xp,yp\Xq,yq)6i92---9N-r-lS%-r(^s+l---SN-l , .n\ 

Triuj yp,Xp\xq,yq) ■ ■ ^ (43) 

T(uj '-Xp,yp\xq,yq)T[Xp,yp\Xq,yq) 

^, , s , f{xp,yp\uj^xq,yq) eX ■ ■ ■ eLi Olds+i ■ ■ ■ Sn-2 

Ts{xp,yp\uj yq,Xq) ^ — ■ ^ (44) 

T{Xp,yp\uj '-Xq,yq)T(Xp,yp\Xq,yq) 

F{Xp, yp\Xq, yq )T{ujXp,yp\ujXq,yq) ^ ^Q^_^iQ^f (45) 

f{Xp, yp\ujXq,yq) f{LUXp, yp\Xq,yq) 

Again, in these relations as written, p,q both lie in the central domain Dq. However, we 
can now analytically continue to any Riemann sheet and use the above automorphisms to 
obtain 

TwVp.q ~ fu-^p,q O162 ■ ■ ■ 6]V-r-l ^AT-r^AT-r + l ' ' ' ^AT-l/ {TjJ-lp^q Tpq) , 

and S are the same as in ]lo|, while R, U occur in . 
Note that in general Up is not obtained from p by staying on the same sheet and simply replacing Xp by u)Xp 
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Tp,u=vq ~ Tp,u=q doOi ■ ■ ■ Og^i ■ ■ ■ O'n-iI {Tp,U-^qTpq) 



(46) 



TpqTup^Uq/{Tp^UqTup,q) ~ (SjV-l/So) 

for all p, q. Note that in these relations r, s must both have values in the set {0, 1, • • . , A^ — 1}. 
Substituting the form ( ^ ) of Tpq, we obtain the exponent relations 

-eu^VpAk + r) = ev^p,q{k + r) - - 1) - epq(k) + [^±ii!llHW-) j _ ^ ^ 

- ep,u^vq(k ~3) = ep,u^q{k - s) - + 1) - e,,(fc) + - [^1 - [#] - (47) 

^pqik) + eup,uq{k) — eup,q{k + 1) — ep^uq(k — 1) = 0po-(i5fc,o — 5fc,]v-i) • 

We have used Mathematica to verify for TV = 2, . . . , 12 that these equations are indeed satisfied 
by (I34I). They are explicitly periodic in r, s, of period A*', so are true for all integers r, s. 
For future reference, we note that the half-exponents for Tp^vq are 

(■p,U''q[k — a) = '4>{p,<J,k,a\m,n) , (48) 

where 

'ilj(p,a,k,a\m,n) = gt+i -~ Qk-i + np^k - cr{mk + nik+i) + (-1)" {riN-a - no - nik-a)- 

aFp-k+<y,c+ pil-cr){Sk,c-l) ■ (49) 
Also, those for and Tu''p,q are 

^u<^p,q{k + Q?) = Qk+i - gk-i — mk+a + p{n-k + "-i-fe) + ( — l)''(mo - niN-a + n-k-a) + 

+ p{l - a){5k,^c. - 1) ■ (50) 



7 The function r2(j9, 

In the derivation of Tpq given in (Q) and ([^), an important role is played by the auxiliary 
function T2(p,q). Noting that the T{xq,yq) in ([^) is Tpq, while T{L0Xq,yq) is Tp^uq, it follows 
that equation (25) of (M]) can be written, for all Riemann sheets, as 



T-2(p,g) 



{Vp - UJXq){tp ~ tq) 



Up {Xp Xq) 



Tr 



p,Uq 



writing T2{tq) in (f^) as T2{p,q). Also, in (53) of ^ we obtain the result 



logr2(p,g) = — 



1 + Ape'" 
1 - XpC'O 



log 



A{e) - UJtq 



Vp 



de 



for p,q both in the central domain ©0, where |Ap|, |Ag| < 1. Here 



m = 



1 - 2k' cos e + k'' 



,2\ 1/iV 



Setting j = Af in (14) of (|j) and using (21) and (23) therein,we find 

N{yp-Xq)(yp-yq){t^ -t^) 



TdqTi 



pq J- p,Vq 



Up {xp 



')iy? -y^)r2{p,U-^q) 



(51) 



(52) 



(53) 



(54) 



for all Riemann sheets. The LHS is the free energy (per double site) of a model with vertical 
rapidities p and alternating horizontal rapidities q, Vq. This is the general "superintegrable" 
model discussed in jl^, so we see that to within simple known scalar factors, l/r2(p, f/~^g) 
is the free energy of this model. Indeed T2{p,U~^ q) itself is the free energy (apart possibly 
from simple scalar factors) of the "inverse" model introduced in ||l^ , while r2 {p, q) is the free 
energy of the model defined in (3.44) - (3.48) of Q (with j = 2 and A; = 0). 
We can take (|5l|) as the definition of T2(p, q). It follows at once that 



T2 (p, q)T2 (p,Uq)--- T2 {p, U q) 



{Vp ^q )(^p iq ) 



(55) 
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Also, considered as a function of tg, T2{p,q) only has a single branch cut, from uj~^ri to 
uj~^ jr]. Across the other A'^ — 1 potential branch cuts it is in fact an analytic function of tq. 
Together with (|5^), this implies that 

T2{v,V^VV-\) = , (56) 

where 

(Going from q to U^VU~^q takes one from Do to the neighbouring domain D^, while leaving 
tq unchanged.) 

One can also verify directly from ( js^ ) that 

T2ip,q)T2{U'VU~'p,q) = {uj~'tp ~ ujtqf /vp . (58) 
Exponent relations 

The exponents of T2{p, q) are simpler than those of the free energy function Tpq, from ( [si] ) and 
(§) , 

r2(p,g)~e,^f,,/Ti,,,, ~ JKefe/e,)''-*"' , (59) 

where 

(^pq(k) = epq{k) — ep^Uqik — I) — 4>pa5k.O 

= (-1)'' (mfc_i - rrifc + no - njv_i - p(5fc_i) . (60) 

We see that the bilnear terms gk do not occur in Ipq (k) , leaving only terms that are fully linear 
in the rrij and rij. 

The equations (^5|), (^6|), ( ^ ) imply that the exponents Ipqik) satisfy 

ipq{k) + ip,Uq{k - 1) H h ep,(7lV-l,(fc — A + 1) = — 4>pa , 

~ip,u-vu--q{k) = epq(fc) + 5»,jv-i , (61) 



for all p, cr, k, mo, . . . , mjv-i, no, . . . , njv-i. Indeed we find, using (|6o| ) and (^, that this is so. 
The half-exponents for Ip^u^q are 

^p.Vqik — a) = ( — l)°'(mfc_„_i — mfc_a + n_a — n-of-i — p5fc,a+i + o-(5ci,jv-i) . (62) 

8 "Sufficiency" of tiie rotation and inversion rela- 
tions 

As is consistent with series expansions, the function Tpq is non-zero and analytic in the central 
(physical) domain Cq. Given this, the recursion relations (^) - ^ are certainly sufficient to 
determine the Riemann surface on which Tpq lives, that Tpq is meromorphic on this surface, 
and to give the orders (exponents) epg(fc) of all its zeros and poles on every sheet. 

This goes a long way towards defining Tpq. To complete the description one needs to 
establish that the ratio of two such functions with the same zeros and poles has some periodicity 
property from sheet to sheet so that it is bounded over the whole surface. It is certainly entire, 
so by Liouville's theorem it would then have to be a constant. Such constants can usually be 
fixed from special cases. 

We shall not discuss this problem of completing the description further herein, but will 
suppose that it can be done. Here our concern is to see if the weaker rotation and inversion 
relations can be used to determine the exponents epq(k). 

More specifically, the rotation and inversion relations, together with the analyticity prop- 
erties in the central domain, are known to be sufficient to determine the free energy for the 
two-dimensional lattice models with the "rapidity difference property" . ^ 



^This property implies that such a model can be parametrized in terms of single-argument Jacobi elliptic 
tions, which fixes the Riemann surface for the free energy. 
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The chiral Potts model does not possess the rapidity difference property, and there is no 
parametrization in terms of single- argument Jacobi elliptic functions. Nevertheless, we can ask 
whether its rotation and inversion relations, together with some simple and plausible ansatz, 
are sufficient to determine the free energy exponents epq{k). This is the question we address 
in the remainder of this paper. 

The relations 

From equations 39, 40 and 10 of [Q, the rotation and inversion relations are 

Tq^Rp = Tpq , TqpTpq = Tpq , (63) 

where 



N{xp~xq)iyp-yq){t^ -t^) 
(x^ -x^){y^ -y^)(tp-tq) 

From these we can deduce that 



Tjip^qTpq — 'i'Rp^q , Tp^RqTpq — ^pq • (6^) 

These two equations can be obtained from (^, by setting r = s = and replacing either Xp 
or Xq by uXp or uXq. They can also be obtained by setting r = s = N — 1 and replacing either 
p or flby Rp or Rq. There is therefore an overlap between ( |63| ) and the recursion relations 
®' ®' latter can not be deduced from the former. 

Using ^ and @, (||) become 

fq,Rp ~ Tpq , Tqpfpq ~ (^1^2 ' ' ' ^JV-l)^ . (66) 

On a given Riemann sheet of p-parity and type p, r, and g-parity and type a, s, specified 
by the integers mo, . . . , m]v-i, no, . . . ,nN-i, let the order of the zero when tp — co'^tq be 
2 (_i)''+°"e(p, CT, fc, r, s|m, n), where m = {mo, . . . , mjv_i} and n = {no, . . . , njv_i}. Then, 
analogously to (p3|), we can write 

Tp, ~ n (^^/^^)''"""'"'"""'"' ' (67) 

the product being over all zeros (and poles), and all Riemann sheeets. 

Are the rotation and inversion relations (pq ) sufficient to fix e{p, a, k, r, s|m, n) as the epq{k) 
given by (B4j)? First we must note that there are severe self-consistency restrictions on the 
exponents e(p, a, fc, r, s|m, n), for any meromorphic function on the Riemann surface. 



Consistency 

Consider some particular zero, at tp = co^tq, on some particular sheet of types r,s. For 
definiteness, take the parities p,a to be zero. Then (since yp ~ ^^,yq — i^") the zero is at 
Xp = u'^+'-'^a;,, yp = w'-'y,. 

Now move p and q to adjacent sheets of types r',s', respectively, so their parities both 
become one. Now consider the zero at tp = tq. This must be when Xp = ui^ ~" Xq and 
yp = uj'' ^" yq. li k — k' — r + r' — s — s' , this is the same zero as the one on the previous 
sheet: we have simply followed it from one sheet to the next. Its exponent (an integer) must 
be the same, so 

e{p,(j,k,r,s\m,n) = e{l ~ p,l — (j,k,r , s \m ,n) , (68) 

provided k = r + r' — s — s' and m', n' are given by (^|), (0)- We have only established this 
condition for p = a — 0, but the same result is obtained for all p, a. 

This is a very strong condition on the exponents e(p, cr, k, r, s\m, n). It must be true for all 
values of p, a, r, s, r' ,s' , and all values of the integers mj,nj. Furthermore , it must be true for 
any meromorphic function on the Riemann surface. Hence it is true for the exponents lpq{k) 
of T2{p,q) as well as those of Tpq and Tpq. 
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Rotation and inversion 

From ( p?! ) - (p2[), exhibiting the dependence of 6^^, 6u on p, q: 

[o.\p = [e-.U ' - [^-A,q ' (69) 

= p-.-i],, - = ■ (70) 

Together with (|6^), (|6^) and (^^, these imply the relations 

e{p,a, k,r, s\m,n) — —e{a,l — p,—k—l,s,r + p\n,Rm) , (71) 

e{p,a,k,r,s\m,n) +e(a,p,~k,s,r\n,m) = (ppa{Sk,o ~ 1) , {T2) 

true for all p, a, k, r, s, m, n. Here Rm is the result of the operation R on the set m, so from 
(H) it follows that {Rm)j = ruj-i + p(l - 5j, 1). 

Other elementary relations 

There are no zeros or poles of Tpq or Tpq in the central domain Do, which is when m = n = 0, 
so 

e(0,0,A;,0,0, |0,0) =0 . (73) 

The model also possesses a reflection symmetry [^. Let S be the automorphism deflned in 
(p6[), (|40[). Then the Boltzmann weights satisfy Ws9,Sp(n-) = Wpq{n), W sq,Sp{n) = 'pypq(— n). 
Replacing p, q by Sg, Sp therefore reflects the lattice about its SW-NE axis. This does not 
change the partition function, so Tsq,Sp = Tpq. 

It also leaves Qij unchanges, while replacing &h by Qji. The product over i,j in ( ps| ) is 
symmetric in i, j, so Gpq is unchanged and from ( pO| ) 

Tsq,Sp = Tpg • (74) 

From (|^ it follows that 

e{l — a,l — p,k, —s,—r\Sn, Sm) = e{p,cj,k,r, s\m,n) , (75) 
where Sm is 5* acting on the set m, so from ( po[ ) {Sm)j = — m]v+i-j. Similarly for S*?!-. 

A bilinear ansatz 

To give the free energy, the rotation and inversion relations must always be supplemented 
by analyticity assumptions, We know from the result ( p^ that e{p,a, k,r, s\m,n) is linear in 
the rrij and rij separately, and hence bilinear in their combination. This seems to be a basic 
property that one may have expected, so we assume that there exist coefficients A, B, C, D 
such that 

JV-l 

e(p, cr, k, r, s\m, n) = A[p, a, k, r, s) + [B{p, o, k, r, s\j)mj + C{p, a, k, r, s\j)nj\ + 

3=0 

JV-l JV-l 

,j)minj (76) 

i = J = 

for all integers rrij , rij . For generality we allow the coefficients to explicitly depend on the 
types r,s of the Riemann sheet. 

There are thus 4A''''(A'^+ 1)^ unknown coefficients to determine. Rather than proceed fully 
algebraically, we have looked at modest values of A'^ (not bigger than 12) and used Mathematica 
(with mj,nj arbitrary) to determine whether the number of variables fixed by the various 
equations. Thus the following remarks are extrapolated conjectures from our observations. 

We take A'^ > 3. As we noted above, the Ising case A'^ = 2 is special in that then the nij 
and rij are not linearly independent, but satisfy mi_p = ni-o- — 0. All the solutions we write 
down are valid for N = 2, but there may be other solutions, so that our remarks concerning 
the number of solutions may not apply. 

First we substituted the ansatz ( |76[ ) into the consistency condition ( |68| ) and observed that 
the number of undetermined coefficients reduced to 2A''^ ( A'' + 1) . Since one can fix A;, r — s and 
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in this condition, this means that for each such case there are just TV + 1 unknown 
coefficients, i.e. N + 1 linearly independent functions satisfying ( jj^ ) . 

We know what these functions are. The functions Ty.U'^a, Tua,, ,, T2(p. U"q) are all mero- 
morphic on the Riemann surface. Their exponents, given 

by ®, (|o|), (H), must therefore 
satisfy (|6^. Taking a = 0, . . . , A'' — 1, we obtain solutions of (S^)- They cannot all be 
linearly independent, but those in the first set certainly are, so that gives us A*' solutions. The 
remaining one is a constant - independent of the rrij and rij. 
It follows that e(p, o, k, r,s\m, n) must be of the form 

JV-l 

e(p, a, k, r, s\m, n) = 7(p, a, k,r — s) + ^ ' c(p, a, k,r — s|ci)i/'(p, a, k, a\m, n) (77) 

where the coefficients 7, c satisfy 

7(p, o", fc,A) — 7(f — p, 1 — CT, fc, fc — A) , c(p, CT, fc,A,a) — c(I — p, 1 — o", fc, fc — A, a) , 

so there are 2N^{N + 1) independent coefficients 7, c, as yet undetermined. The consistency 
condition (^) is now satisfied. 

For N odd, the rotation relation (|7l|) reduces the number of undetermined coefficients to 
N{N + 1). Then the inversion relation ( JT^ ) fixes all the coefiicients, giving the solution (^). 
The extra relations ([ts]), ( [ts] ) are then satisfied. 

For A'^ even and greater than 2, (^l|) reduces the number of undetermined coefficients to 
N{5N + 6)/4. Then ^) further reduces the number to [(iV + 4)/4] (so for iV = 4, 6, 8, 10, 
12 the numbers are 2, 2, 3, 3, 4). The other relations ([73[), ( [tH} ) are satisfied by this solution, 
so do not reduce the number of undetermined coefficients any further. The solution is then of 
the form 

epg(fc) + h{p — a - k + 2r — 2s)/(p, a, k\m,n) , 



where e-pq{k) is given by ( |34[ ) and 

JV-l 

f{p,a,k\m,n) = 5fc+i - 5fe-i - mfe+a + rip-fc + ^ [(-l)°'~'''~'^ma - (-l)"+*~''na] . 

These functions have the weak difference property that they are unchanged (for A'^ even) by 
incrementing all of mo, . . . ,mjv-i,no, . . . ,njv-i by the same arbitrary integer, so no further 
restrictions can be obtained by making this requirement. 
The coefficients h{j) are integers, subject only to the constraints 

h{j) = h{N - j) = h{N + j) , and h{j) = if j is odd . 

It follows that just [{N + 4)/4] of them (all with j even) remain undetermined. 

To summarize: if A' is odd, the rotation and inversion relations, together with the consis- 
tency condition (^) and the bilinear ansatz (^^, are sufficient to determine the exponents 
e{p,(j,k,r,s\m,n). Surprisingly, one does not need the analyticity and non-zeroedness of Tpq 
in the central physical domain 7?o, i.e. the relation (|73|). 

The same is true for N even and p — a — k odd. In both these cases we find that there is 
no explicit dependence of e{p, a, k, r, s|m, n) on the types r, s of the Riemann sheet. 

For A' even and p — a — k even they are not quite sufficient, but they leave only [(A'^ + 
4)/4] parameters to be determined. If one also (guided by the other cases) assumes that 
e{p,a,k,r,s\m,n) is not explicitly dependent on r or s, then there is only one undetermined 
coefficient left. 

As we remarked above, a knowledge of the orders of the poles and zeros of a meromorphic 
function does not by itself fix the function, but it goes a long way towards it. For instance, 
one might observe that the orders tpq{k) satisfy the relations ( p7| ) and then guess the full set 
of recursion relations (g) - (||), and hence (|) - (0). One could then test these relations from 
series expansions. They imply the vital Assumption 2 of so one could then use the method 
of Q to obtain the result (67) therein. 
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9 Summary 



We have shown that the partition function per site T^q lives on a Riemann surface with an 
infinite number of sheets, each sheet corresponding to a point an (2A'' — l)-dimensional lattice. 
However, one should note that adjacent sheets need not correspond to any simple geometric 
definition of adjacent points: the adjacency rules are given by (1^), (p^). 

It is a meromorphic function on this surface, with zeros and poles only when tp = uj^tq 
for = 0, . . . , A'^ — 1. The orders of these zeros (the negative of the order of the poles) are 
bilinear in the integers mo, . . . , n]v-i that specify the sheet. If we assume this bilinearity 
(with coefficients that may depend on k and the parities and types of the Riemann sheet), 
then for A'^ odd they can be obtained from the rotation and inversion relations. For A'^ even 
this procedure does not uniquely fix the orders, but does determine them to within a small 
number [(A'^ + 4)/4] of free parameters. 

The integers mo, . . . ,njv_i (more precisely the differences of the ruj, and of the Uj) are 
connected with the variables of the hyperelliptic paarametrization of the chiral Potts model 

(with u,v therein interchanged). We hope to discuss this point in a later paper, and at 
least present an elliptic function expression for T2(tq) in the case A'^ = 3. 

A significant motivation for this work has been the still outstanding problem of obtaining 
the spontaneous magnetization (order parameter) of the chiral Potts model. There is an 
elegant conjecture (eqn 3.13 of |^| , eqn. 1.20 of |Q| , eqn. 15 of |^|) for this property, which 
is almost certainly true, but has not been proved. Following the method of Jimbo, Miwa and 
Nak ayas hiki |l6| , the author has derived functional relations for a generalized order parameter 
[0, jl^. These have a similar structure to the inversion/rotation relations for the free energy. 
If one could solve them, then one would have verified the conjecture. 

It was therefore the author's hope that the techniques of this paper could be applied 
to solving the order parameter relations. It must be admitted that preliminary results are 
not encouraging. If we assume that the function lives on the same Riemann surface as Tpq 
(wwhich is not obvious), and that the orders (exponents) of the poles and zeros when — tq 
are bilinear in the mj,nj, then for A'^ odd we find the solution is unique, but is merely the 
"wrong solution" we obtained in equation 72 of |l^. So we are no further forward. 

It is possible that this function has poles and zeros other than those when tp = tq : we 
have found some suggestion of this in preliminary low-temperature expansion calculations, 
and hope to discuss this in a later paper. If so, then the order parameter function may be 
much more complicated than Tpq, and the simple ideas we used here may need considerable 
expansion. 
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